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ABSTRACT 

We present a spectroscopic and photometric analysis of the remarkable mas- 
sive binary system, BD+36°4063. We argue that the visible ON star is under- 
going a rapid mass transfer episode that results in a thick and opaque disk that 
surrounds and renders invisible its massive companion. A comparison of the 
projected rotational velocity and the orbital semiamplitude of the visible star 
indicates a mass ratio near unity. Models for conservative mass transfer show 
that the equal mass state occurs at the point of minimum separation, and we 
argue that BD+36°4063 may represent the first system identified at this rapid 
and rare stage of massive binary evolution. 

Subject headings: binaries: spectroscopic — stars: early-type — stars: evolution 
- stars: individual (BD+36°4063) 



1. Introduction 

Massive close binaries are destined to un dergo extreme trans formations once they begin 



to exchange mass and angular momentum (ILanger et al.ll2008l ). If mass transfer is con 



servative, then the system separation decreases as mass transfer progresses, and the orbit 
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dimensions reach a minimum once a mass ratio of unity is reached. The syst em then experi- 



ences a critical phase where the mass transfer rates may exceed 10" 4 M Q y" 1 (jWellstein et al.1 



20011 ) . Subsequently, the orbit will widen after the mass donor becomes the lower mass ob- 
ject, and a slower and more extended mass transfer stage will begin. The gainer accepts 
both mass and angular momentum, and it may quickly reach a state of critical rotation 
where gravitational and centripetal accelerations balance at the equator. This may limit 
further mass accretion and l ead to the formation of a thic k, opaque, and mass losing disk 
surrounding the gainer star (INazarenko fc Glazunovall2006l ). 



This stage o f binary evo l ution is prob ably represented by the W Serpentis class of 
massive binaries (lPlavedll980l ; iTarasovl |2000| ) . These binaries contain a visible, low gravity, 
donor star that fills its Roche lobe, while the higher mass, gainer star is invisible in the 
composite spectrum because it is immersed in a thick disk. The first direct image of such 
a disk around a gainer star was recently made with the CHARA Array for the best known 
member of the class (3 Lyrae (IZhao et al.ll2008l ). The shortest period W Ser binary known is 
RY Scuti (P = ll. 



( iGrundstrom et al 



d). I t consists of a 7M , 09.7 Ibpe star and a 28M© hidden companion 
200 7|), an d the entire system is surrounded by a 2000 AU wide gas and 



dust torus (ISmith et al.ll2002l ). The known W Ser systems are probably experiencing the slow 
mass transfer stage after mass ratio reversal, and no W Ser system has yet been identified 
at the close, rapid mass transfer stage. 

The mi ssing link may be filled by the subject of this paper, BD+36°406 3. The star 
(ON9 .7 lab; IWalborn fc Howarth 120001) is located in the Cyg OBI association ([Humphreys 



1978 ). It was noted by Mathysl ( 1989 ) as a key example of an ON star, showing spec- 
tral featur es indicating the N-enrichm ent and C-depletion characteristic of CNO-processed 
elements (IWalborn fc Howarth! 120001 ) . Many ON stars are members of binary systems 
( Bolton fc Rogers 19781 ) where mass transfer may remove the H envelope and reveal CNO- 
processed gas in the atmosphere. Thus, i t is very pertinent th at BD+36°4063 was discovered 
by Howarth (IWalborn fc Howarth! |2000| ; lHarries et al.l |2002| ) to be a single-lined, spectro- 
scopic binary with a period of 4.8 d. In a recent meeting contribution^], Howarth presents 
radial velocity and light curves, and he shows that although the companion is massive, its 
absorption features are completely absent from the observed spectrum. 

We obtained blue spectra of this and several other faint O-stars in a search for binaries, 
unaware at the time of Howarth's detection of binary motion. Here we confirm the orbital 
period and present preliminary orbital elements for the binary. We argue that the visible 
star is filling its Roche lobe and that the high rate of mass transfer results in a dense disk 
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that blocks a direct view of the massive companion. We conclude with a brief discussion of 
system parameters and evolutionary status. 



2. Observations and Orbital Elements 

Seven spectra of BD+36°4063 were obtained with the Kitt Peak National Observatory 
(KPNO) 2.1 m telescope from 2008 November 15 to 21. These observations made use of the 
Goldcam spectrograph with grating G47 (831 grooves mm -1 ) in second order with a CuSO.4 
order sorting filter. The detector was the T3KC CCD (a 3072 x 1024 pixel array with 
15 x 15 yum pixels), and the resulting spectra have a resolving power of R = A/AA = 2400 
as measured from the HeNeAr comparison lines. Exposures were usually 600 s in duration, 
leading to spectra with a S/N = 200 per pixel in the continuum. The wavelength range is 3942 
- 5032 A with a wavelength calibration accuracy of ~5 km s _1 based on the rms scatter of fits 
to comparison lines and on the variance of multi-night measurements of the spectra of other 
stars. The spectra were extracted and calibrated using standard routines in and then 

each continuum-rectified spectrum was transformed to a common heliocentric wavelength 
grid in log A increments. 

Only one set of spectral lines was readily apparent in these spectra, so we deter- 
mined radial velocities by cross-correlati ng each spectrum with a model spectrum from the 



TLUSTY/SYNSPEC BSTAR2006 grid (lLanz & Hubenvi l2007h . We selected a model for 



T e ff = 28 kK, log g = 3.0, and mic roturbulent veloci t y of 10 km s \ par ameters that are 



typical for the star's classification ( IRyans et al.ll2002l ; Martins et al.ll2005l ). although small 



variations in these stellar parameters do not significantly alter the resulting velocities. We 
selected a N-enriched "CN" model to match better the line spectrum. The model spectrum 
was transformed to the observed grid by integration and was then convolved with functions 
to account for rotational broadening (§3) and instrumental broadening. A number of spectral 
regions that contain ISM features or emission lines were omitted from the cross-correlation 
sample (all hydrogen lines, the diffuse interstellar band at A4428, and He II A4686). Table 1 
lists the date of observation, spectroscopic phase, radial velocity and error, and the observed 
minus calculated residual from the fit (below). 

Our seven day time span is much too short to derive an accurate orbital period, so we 
estimated the period from the ellipsoidal variations in flux related to the tidal distortion 
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of the star. Two sets of photometric measurements exist from all sky survey experiments. 
The first set of 111 Cousins In measurements were made between 2003 and 2007 by The 
Amateur Sky Survey (TASS9; broege et alfaod ). The second set of 84 points were made in 



Amateur bky Survey { lAbqj lUroege et ai.ll/UUbl ). ine s econd set ot o4 points were made m 
1999 with the Northern Sky Variability Survey (NSV^I IvVozniak et alil2004h . Because the 



NSVS measurements record a broad spectral range (4500 - 10000 A), we simply subtracted 
a constant value of 1.099 mag to bring their mean into coincidence with the mean of the 
/,7-band results. Despite differences in wavelength intervals, both the amplitude of variation 
and the ephemeri determined from each data set matched within uncertainties. A discrete 
Fourier transform of this combined time series immediately showed evidence of half the 
orbital period. A double sinusoid light curve is exhibited in one orbit by tidally distorted 
stars, and for our data this yielded a period of P = 4.8126 ± 0.0004 d and an epoch of ON 
star superior conjunction at T^c = HJD 2,452,448.60 ± 0.08. 

We then determined the remaining ele ments using a fit of th e veloc ities made with the 



non- linear, least-squares fitting program of iMorbey fc Brosterhud ( 119741 ) . The weights were 



set to unity for most measurements (since they have comparable measurement errors), but we 
increased the weight of the single datum in the negative branch to four (to balance the four 
measurements in the positive branch) and we decreased the weight of the two conjunction 
phase measurements to 0.5. We fixed the period at the photometrically derived value and 
then solved for the circular orbital elements: systemic velocity 7 = —17 ±3 km s _1 , velocity 
semiamplitude K = 163 ±3 km s -1 , and epoch of superior conjunction of the ON star Tsc = 
HJD 2,454,787.26 ±0.03. Trial elliptical solutions did not improve the fit. The rms of the fit, 
6.6 km s _1 , is larger than the formal errors but comparable to what we find in observations 
of other O-stars from this run. The spectroscopic epoch occurs (0.3 ±0.2) d earlier than the 
prediction from the light curve ephemeris, which may indicate that the period is decreasing. 
The possibility of a measurable period change offers us an important diagnostic of the mass 
transfer rate that must be confirmed in future observations. 



3. Binary Properties 

The most surprising feature of the orbital solution is that the mass function is large, 
f(m) = (2.18 ± 0.12)M Q , suggesting that the companion is a massive star. However, there 
is no clear evidence of absorption lines of the companion in the individual spectra. The 
expected radial velocities for any companion lines depend on the assumed mass ratio. We 
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can estimate the mass ratio by considering the projected rotational velocity V sin i of the ON 
star. If we assume that the ON star fills its Roche lobe and rotates synchronously, then the 
ratio of V sin 2 to semiamplitude K is a monotonically i ncreasing function o f Q = Mi/M 2 
(mass of the ON star divided by that of the hidden star; iGies &: Bolton! 1 19861 ) 



I^sini 
K 



(Q + 1)$(Q) 



where is the fractional Roche radius of the ON star (jEggletonlll983l ). Thus, we can use 
a measurement of Vsini to determine the mass ratio. The resolution of our spectra is just 
adequate for this task. We made measurements of the FWHM of the deep and relatively 
unblended line of Si IV A4088, whose profile should be dominated by rotational broadening. 
The mean width of the line is FWHM = 2.95 ± 0.11 A, which is significantly larger than 
the instrumental broadening measured in comparison lines near this wavelength, FWHM 
= 1.90 ± 0.02 A. We then created synthetic spectra for a grid of tes t values of V sin % by 
convolving the model spectrum with a rot ational broadening func tion (lGrayll2005l ) for a lin- 
ear limb darkening coefficient of e = 0.37 (jWade fc Rucinskilll985l ). It should be noted that 
varying the linear limb darkening coefficient even by a factor of two has no effect, within 
uncertainties, on V sin i measurements. These model profiles match the observed FWHM for 
Vsini = 126 ±15 km s~ x . Although it is possible that some of the apparent line broadening 
is due to macroturbulence in the atmosphere, a test using a presu med macroturbule nt broad- 
ening with £ = 50 km s _1 (at the high end for similar supergiants; iRyans et al.ll2002l ) resulted 
in the same V^sinz because the instrumental broadening is so much larger than the expected 
macroturbulent broadening. Then the relation above leads to a mass ratio Q = 1.02 ± 0.17, 
i.e., the stars are about equal in mass. 

We see no evidence of absorption lines moving in anti-phase with a velocity comparable 
to that of the ON star. However, Figure 1 illustrates the orbital variations of the H/? line 
in which an emission component does appear to share the orbital motion of the companion. 
We made a preliminar y Doppler tomographic reconstruction of the two spectral components 
( IBagnuolo et al.lll994j ) assuming the mass ratio given above, and we found evidence of similar, 
anti-phase moving emission lines in H7, He I AA4387, 4471, 4713, 4921, and He II A4686. We 
suggest that these emission features form in a disk surrounding the gainer star. Significant 
emission lines also are observed at Ha (see Howarth's conference poster) and in the near- 
IR at He II 2.058 fim and Br7 2.166 /an jHanson et all Il996l : iTamblvn et~aD Il996h . A 
comparison of the line depths in the reconstructed spectrum of the ON star with those in 
the model indicates that the continuum flux from the disk is faint, F 2 /Fi » 0.1. Note that 
the absorption part of H/? appears deeper at conjunctions. The same deepening was observed 
in other lines that strengthen in slightly cooler atmospheres (Si III, O II). Similar changes 
are found in RY Scuti and are probably related to the tidal extension and gravity darkening 
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of the donor star (IGrundstrom et al.l 120071 ). 



The example of RY Scuti shows that mass transfer to a thick disk can lead to sys- 
temic mass loss from th e binary and the formation of circumbinary gas and dust struc- 
tures (ISmith et al.l 120021 ) . There is some indication that an infrared excess from such cir- 
cumbinary material is present in the spectral energy distribution (SED) of BD+3 6°4063. 



We show in Figure 2 th e available fl ux measurements based u pon Johnson UBV (IHiltner 



1956 
2006 



Colina et al. 


1996^ 


Bessell et al. 


1998* 



1998|), and 2MASS JHK S magnitudes flCohen et alJl2003 



1993) . Stromgren JCrawford [l~9~75l: lGravl[l9~9~8l) TASS V. In. JProeee et al 



Skrutskie et al. 



2006), plus mid-IR fluxes from the Spitzer/IRS post basic calibrated data archivqj ( IHouck et al 
2004). We fit these ob served fluxes with the BSTAR2006 flux model to determine the redden- 
ing (jFitzpatricklll999l ) E(B — V) — 1.28 ± 0.06 mag and ratio of total-to-selective exti nction 
Ry = 3.50 ±0.10 (consistent with previous reddening estimates; iPatriarchi et al.ll2003l ). The 
limb darkened, angular diameter of the ON star is Old = 97 ± 8 //as (after correction for a 
companion flux contribution of F2/F1 = 0.1). The resulting radius - distance relationship is 
Ri/Rq = (10.4 ± 0.8) d(kpc). There does appear to be a flux excess in the SED at wave- 
lengths > 8/j,m that may result from systemic mass loss. Could stellar winds account for this 
excess? Because the BSTAR2006 grid does not take winds into account, we also fit an SED 
from a CMFGEN model atmosphere (IHillier &: Miller! Il998l ) with similar stellar parameters 
and a wind. The absolute fluxes from the CMFGEN model were systematically 5% higher 
than those from the BSTAR2006 SED over the optical to near-IR range. However, there is 
no evidence of a marked IR excess near 10 /im from the wind in the CMFGEN model SED. 
This suggests that the observed excess may result instead from binary mass loss. 

We summarize the available constraints on the masses of the stars in Figure 3. The mass 
ratio range derived from Vsini/K is indicated by the lines of constant slope. Next, we can 
use the observed light curve to constrain the orbital inclination. The ON star probably fills 
out its Roche surface, and consequently the photometric variations are due to the star's tidal 
distortion (geometric shape and gravity darkening). The amplitude of the ellipsoidal light 
curve depends on the degree of Roche filling (assumed complete) and the orbital inclination 
(larger at higher inclination). Ideally we would create a model that includes bo th a Roche 
distorted mass donor and a disk surrounding the gainer (jDjurasevic et al.ll2008l ). but given 
the small flux contribution of the disk and the relatively large errors in the photometric 
data, we made an approximate model assuming that the companion is a small spherical 
object (with the same temperature as the ON star) that acts only as a mild flux dilution 
source in the light curve of the tidally distorted ON star. We used the GENSYN binary 
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code (IMochnacki &: Doughtyl Il972l ) to create model Ic light curves for a grid of orbital 
inclinations, assuming that the ON star fills its Roche surface and that the companion has a 
radius that yields a monochromatic flux ratio, F%/F\ = 0.1 (matching that of the preliminary 
tomographic reconstruction). We obtained a best fit light curve (full amplitude ~ 0.13 mag) 
for an inclination % = 49° ±8°. The derived masses based on the spectroscopic mass function 
and this range in inclination are sketched as curved lines in Figure 3. 

We can combine the radius - distance relation derived from the SED with the projected 
rotational velocity (assuming a Roche filling, synchronously rotating, ON star) to form a 
distance - inclination relation, sinz = (1.15 ± 0.16)/d(kpc). The inclination range shown 
in Figure 3 corresponds to distance range of 1.37 kpc (i = 57°) to 1.76 kpc (i = 41°), 
which agrees well wit h distance estimates for the Cyg OBI association, d = 1.25 — 1.83 kpc 
( lUyaniker et al.ll200ll ). We can also associate a stellar luminosity with any specific distance 
(or inclination) from the radius - distance relation. In the absence of mass transfer, we would 
expect the lu minosity to agree w ith the mass - luminosity relation for single stars of its 
temperature (ISchaller et al.lll992l ). so that one position along the constant inclination locus 
would be preferred. However, stars in binaries that suffer mass loss may appear overluminous 
for their mass, so in practice there is an upper limit along the inclination track where the 
star becomes more massive than expected for the lum inosity. These terminal points are 
indicated in Figure 3. Because the ISchaller et al.l (119921 ) models to not account for r otation 
and because rapidly rotating stars may appear more luminous ( lEkstrom et al.ll2008l ). these 
end-point values should be taken as upper limits. The gray shaded region shows the preferred 
mass ranges that fulfill all the constraints, and the optimal fit occurs for Mi ~ M 2 ~ 21M Q . 

The tentative picture that emerges from the spectroscopic analysis is that the ON star 
is transferring gas to a relatively faint disk that surrounds and obscures the companion star. 
The ON star must be very close to filling its Roche surface since we observe evidence of 
its tidal distortion in the light curve and in the deepening of lines at conjunctions. The 
emission features that follow the orbital radial velocity characteristics of the hidden mass 
gainer probably form in the thick disk. Thus, BD+36°4063 shares many features in common 
with the W Ser class of interacting binaries, but since the orbital period is much shorter 
and the mass ratio closer to unity than that found in other W Ser systems, we suggest that 
BD+36°4063 represents an earlier and faster mass transfer stage of evolutio n. If the system 



is coe val with the stars of Cyg OBI (with an age of approximately 7.5 Myr; lUyaniker et al. 



2001), then the mass - radius - age relations for single star evolutionary tracks (ISchaller et al. 



19921 ) suggest that the ON star probably began life with a mass of 25M or less. Thus, the 
ON star has probably lost only a modest fraction of its original mass so far. BD+36°4063 
offers us an important glimpse of binary evolution at its most intense stage, and we encourage 
new observational efforts to probe the system and its environment. 
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Table 1. Radial Velocity Measurements 



HJD 


Orbital 




AV r 


O-C 


(-2,400,000) 


Phase 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


54785.624 


0.660 


126.7 


2.4 


5.1 


54786.576 


0.858 


109.7 


2.3 


-0.8 


54787.557 


0.062 


-61.9 


1.9 


16.7 


54788.558 


0.270 


-180.0 


2.2 


-1.2 


54789.567 


0.480 


-36.7 


1.8 


0.7 


54790.556 


0.685 


132.1 


2.1 


-1.3 


54791.559 


0.894 


77.9 


2.1 


-6.8 
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Fig. 1. — The variations of H/3 A4861 as a function of orbital phase and velocity in the frame 
of the ON star. The intensity between observed spectra in the gray-scale image is calculated 
by a linear interpolation between the closest observed phases (shown by arrows along the 
right axis). The backwards S-curve in the gray-scale image shows the relative velocity curve 
of the companion star. An emission component (bright in the gray-scale image) appears to 
follow the velocity curve for the companion star. 
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Fig. 2. — The spectral energy distribution for BD+36°4063. The plus signs indicate the 
observed fluxes while the solid line shows a reddened version of a BSTAR2006 flux model 
for the ON star. 
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Fig. 3. — A mass plane diagram showing observational constraints on the masses of the 
ON star Mi and the unseen companion Mi- The dashed line shows the mass ratio derived 
from Vsini/K and the adjacent dot-dashed lines give the associated ±lcr error range. The 
solid line shows the relation from the spectroscopic mass function for an orbital inclination 
of % — 49° derived from the light curve, while the dotted lines show the same for % = 41° 
(above) and % = 57° (below). These three lines are terminated where the ON star would 
be underluminous for its mass. The shaded region indicates the ranges meeting all the 
constraints. 



